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Abstract 

£NJ ■ In this article, we explore the ability of direct and indirect dark matter experiments to not only 

detect neutralino dark matter, but to constrain and measure the parameters of supersymmetry. In par- 
ticular, we explore the relationship between the phenomenological quantities relevant to dark matter 
t^. \ experiments, such as the neutralino annihilation and elastic scattering cross sections, and the under- 

lying characteristics of the supersymmetric model, such as the values of /x (and the composition of 
' the lightest neutralino), uia and tan (3. We explore a broad range of supersymmetric models and then 

focus on a smaller set of benchmark models. We find that by combining astrophysical observations 
with collider measurements, fi can often be constrained far more tightly than it can be from LHC data 
alone. In models in the ^4-funnel region of parameter space, we find that dark matter experiments can 
potentially determine to roughly ±100 GeV, even when heavy neutral MSSM Higgs bosons (A, 
Hi) cannot be observed at the LHC. The information provided by astrophysical experiments is often 
highly complementary to the information most easily ascertained at colliders. 
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I. INTRODUCTION 



A great deal of effort has been directed to developing methods of detecting particle dark 
matter. Over the past years and decades, numerous studies have been conducted to assess 
the prospects for these various techniques 0. If dark matter consists of neutralinos, or an- 
other weakly interacting particle with a TeV-scale mass, it is likely that one or more of these 
techniques will make the first detection of dark matter particles within the next several years. 

Of all of the candidates for dark matter that have been proposed, none has received as much 
attention as the lightest neutralino in models of supersymmetry. Supersymmetry is theoretically 
attractive for a variety of reasons. Among the most compelling is its ability to provide a natural 
solution to the hierarchy problem , and a common scale for the unification of the forces of the 
Standard Model From the standpoint of providing a dark matter candidate, the lightest 
neutralino is naturally stable by virtue of R-parity conservation pj, and in many models is 
produced in the early universe in a quantity similar to the measured density of cold dark 
matter 0. 

Astrophysical techniques for detecting neutralinos include direct and indirect detection ex- 
periments. Direct detection experiments attempt to observe neutralinos scattering elastically 
off of target nuclei. Indirect detection experiments, in contrast, attempt to detect the anni- 
hilation products of neutralinos, including gamma-rays, neutrinos, positrons, anti-protons and 
anti-deuterons. 

Direct and indirect detection measurements are of critical importance in determining the 
identity of dark matter. Even if collider experiments were to observe a long-lived, weakly inter- 
acting, massive particle that appears to be a suitable dark matter candidate, such experiments 
will never be able to determine whether a particle is stable over cosmological timescales. To 
determine whether the dark matter of our universe is made up of such a particle (either entirely, 
or in part), direct and indirect detection experiments will be needed. 

But looking beyond the mere detection of dark matter, what will these astrophysical ob- 
servations reveal to us about the nature of particle dark matter? In particular, is it possible 
to determine the properties of the lightest neutralino and the corresponding supersymmetry 
model by direct or indirect dark matter detection experiments? In this article, we attempt to 
address these questions. While extracting such information from direct and indirect dark mat- 
ter detection experiments is challenging, it may be possible in many scenarios. We have found 
that direct detection measurements of the neutralino's spin-independent elastic scattering cross 
section, rates in neutrino telescopes, and the brightness of gamma-ray lines (from XiXi ~^ 77 
and XiXi ~^ 1%) are among the most useful astrophysical probes for determining the proper- 
ties of supersymmetry. In most cases, the value of \x (or alternatively, the composition of the 
lightest neutralino) can be constrained far more tightly if astrophysical data is included than 
it can be by LHC data alone. Furthermore, in some models (those in the A-funnel region of 
supersymmetric parameter space) the mass of the CP-odd Higgs boson can be determined by 
astrophysical experiments to within roughly ±100 GeV, even if it cannot be observed at the 
LHC. 



II. THE COMPOSITION OF THE LIGHTEST NEUTRALINO 

In the Minimal Supersymmetric Standard Model (MSSM), the neutral electroweak gauginos 
(B,W) and higgsinos (Hi,H 2 ) have the same quantum numbers and, therefore, mix into four 
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mass eigenstates called neutralinos. The neutralino mass matrix in the B-W-Hi-H 2 basis is 
given by 
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where Mi, M 2 and /i are the bino and wino masses, and the higgsino mass parameter, respec- 
tively. 9w is the Weinberg angle, and tan/? is the ratio of the vacuum expectation values of 
the up and down Higgs doublets. This matrix can be diagonalized into mass eigenstates by 



Mo = N^M Y oN. 



(2) 



In terms of the elements of the matrix, N, the lightest neutralino (x?, or simply x°) is the 
following mixture of gauginos and higgsinos: 



N n B + N 12 W 3 + N 13 Ht + N U H 2 



(3) 



The absolute square of these coefficients are known as the bino, wino and higgsino fractions of 
the lightest neutralino, respectively. The composition of the lightest neutralino is a function of 
four supersymmetric parameters, Mi, M 2 , Li and tan/?. This becomes further simplified if the 
gaugino masses are assumed to evolve to a single value at the scale of Grand Unification, as 
they often are. Requiring that Mi, M 2 and M 3 (M 3 being the gluino mass) evolve to the same 
value at the GUT scale yields the following ratios at the electroweak scale: 



Mi = - tan 2 6 W M 2 « 0.5M 2 , 
3 



M, 



as 

O-EM 



sin 2 9 W M 2 w 3.7M 2 . 



(4) 
(5) 



Adopting this relationship between Mi and M 2 , we plot the composition of the lightest neu- 
tralino as a function of M 2 and \i in figure^for the case of tan (3 = 10. In each frame, the solid 
lines represent contours of constant bino, wino or higgsino fraction, as labeled. The dashed lines 
represent contours of constant neutralino mass, in GeV. For M 2 /2 li the lightest neutralino 
is almost entirely bino-like, whereas for M 2 /2 ^> \x it is a nearly pure-higgsino, containing large 
quantities of both Hi and H 2 . As a result of Mi always being significantly smaller than M 2 , 
the lightest neutralino is never primarily wino-like. This, of course, is not always true if the 
gaugino masses do not obey the GUT relationship of Eq. HJ 



III. NEUTRALINO ANNIHILATIONS 

The annihilation cross section of the lightest neutralino is important for determining its 
thermal abundance in the universe, and for determining the fluxes of neutralino annihilation 
products that may be observed in indirect detection experiments. Furthermore, by studying 
future signals from indirect detection experiments, such as the spectrum of gamma-rays gen- 
erated in the annihilations of neutralino dark matter, the annihilation modes of the lightest 
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FIG. 1: Contours of constant bino, wino and higgsino fractions of the lightest neutralino in the M2-H 
plane. The GUT relation, M\ = | tan 2 OwM^, has been adopted, and tan (3 = 10 has been used. Also 
shown as dashed lines are contours of constant neutralino mass, in GeV. 



neutralino could (in principle) be identified. Neutralino annihilations which generate these sig- 
nals take place in situations in which the neutralinos are traveling at velocities far below the 
speed of light. Therefore, we can focus here on the annihilation modes preferred by neutralinos 
in the low velocity limit. 

Neutralino annihilations can produce a large variety of final states, including pairs of fermions 
and pairs of gauge or Higgs bosons. We will briefly summarize the most relevant aspects of 
these processes here. For a more complete discussion of neutralino annihilations, see Ref.0] 
and references therein. 
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A. Neutralino Annihilations to Fermions 



Neutralinos can annihilate to fermion pairs through three classes of tree level Feynman 
diagrams: s-channel Higgs boson (Hi, H 2 or A) exchange, t-channel sfermion exchange, or 
s-channel Z boson exchange. In the low velocity limit, the exchange of A, sfermions and Z 
each contribute. The cross section resulting from Higgs exchange is given by: 

AG\c f m 2 x0 m 2 w m 2 f tan 2 ft C\ 0A (l - mj/m 2 ^) 1 / 2 

cr Y o y o^ A _+ f fV (v — > 0) = > 1 — — — 5 1 1 — 5— r 1 , (6) 

VjV J ' vr[(4m 2 - m 2 A ) 2 + m 2 A T 2 A ] 



where the sum is over fermion species, Cf is a color factor (3 for quarks, 1 for leptons), m A is 
the CP-odd Higgs mass, and Ya is its width. The expression shown is valid in the case of down- 
type fermions. For up-type fermions in the final state, the tan 2 ft should be replaced by cot 2 ft. 
The neutralino-neutralino-Higgs coupling is given by, C x o x o A (g 2 Ni2 — giNu)(Nu cos ft — 
N 13 smft). This contribution to the annihilation cross section is most large in the case of a 
mixed gaugino-higgsino with large tan ft and with a mass near the A-pole (m^ ~ 2m x o), known 
as the funnel region. 

The contribution from sfermion exchange is given by: 

_ AG 2 c f m 2 m 4 w C 2 Jl - m 2 lm 2 ) l l 2 

a x\x^I^ff V ^^^-2^ n[m 2^ {1 _ m J /m 2 ?) + m 2]2 > 

where, in the small Ni 2 limit, Cj ~ (m/A r 11 A^* 3 taxiOw l^rnw cos ft) + 
3 j\N 13 \ 2 /2m 2 v m x o cos 2 ft)+uif tan 2 9 W |A r n| 2 (0.25+2e 2 +e/) for down-type fermions, and Cj 



(-m f N n N* 4 tan9 w /2m w sin ft) + (m 3 f \N u \ 2 /2m 2 v m x o sin 2 /?)+m / tan 2 6» V y|A r ii| 2 (0.25+2e 2 ,-e / ) 
for up- type fermions. Cf is the charge of the fermion. Unlike in the case of A exchange, this 
process can contribute significantly for either a mixed bino-higgsino or a pure bino. 
The contribution from Z boson exchange is given by: 

f? fC ^(iJv u PH^mi- m }w 

/ 



which scales simply as the square of the difference of the two higgsino fractions, (|A^ 13 | 2 — lA^^) 2 ) 2 . 

Before moving on to neutralino annihilations to bosonic final states, a few comments are in 
order. Firstly, all annihilation channels to fermions strongly prefer heavy fermions in the low 
velocity limit (<jv oc in 2 ). Secondly, interference terms between each of these three channels 
can be important and, in some cases, dominate (see Ref. 0| for the complete cross sections, 
including interference terms). 



B. Neutralino Annihilations to Gauge Bosons 

Neutralinos can annihilate to W + W~ though t-channel chargino exchange, as well as s- 
channel exchange of a Z boson or CP-even Higgs boson. The diagrams involving Z and Higgs 
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bosons do not contribute in the low velocity limit, however. The chargino exchange diagrams 
result in an annihilation cross section given by: 



<J x o x o^ w+w -v{v -> 0) 



±G 2 F m 7 w ((m 2 /m 2 w ) - !) 3/2 
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(9) 



where <p + and 0_ are the angles appearing in the chargino mixing matrices. 

Neutralino annihilations to ZZ can take place through t-channel neutralino exchange and 
s-channel Higgs exchange. Again, the s-channel processes do not contribute in the low velocity 
limit. The neutralino exchange diagrams yield: 



a xWl^zz 
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(10) 



Note that neither of these expressions contain the bino-content of the neutralino. Both channels 
to W + W~ and ZZ are most important for neutralinos with a large higgsino component (or a 
large wino component in the case of iy + iy _ ). 



C. Neutralino Annihilations to Final States Including Higgs Bosons 

Neutralinos can also annihilate to pairs of Higgs bosons, or to a Higgs boson along with a 
W ± or Z boson. Possible final states include H X H X , H X E 2 , H 2 H 2 , AA, AH U AH 2 , H + H~, 
ZHi, ZH 2 , ZA and W ± H T , although processes to H^, H X H 2 , H 2 H 2 , AA, ZA and H+H~ do 
not contribute in the low velocity limit. Of the contributing final states, those including a Higgs 
and a gauge boson (ZHi, ZH 2 and W ± H T ) are generally more important than Higgs-only final 
states, so we focus on those here. 

Neutralino annihilations to ZH\ and ZH 2 can proceed through the t-channel exchange of 
a neutralino, or the s-channel exchange of Z or A, all of which contribute in the low velocity 
limit. Similarly, neutralinos can annihilate to W ± H T through t-channel chargino exchange, or 
the s-channel exchange of H 2 , Hi or A, of which chargino and A diagrams contribute in the 
low velocity limit. 

The cross sections to these final states are larger for higgsino-like neutralinos, but decrease 
more slowly with increasing bino-fraction than annihilations to ZZ or W + W~. In particular, 
annihilations through an A to ZH\ and ZH 2 are proportional to 1 7Vi3 1 2 rather than some 
combination of iVi3 and A?i4 to the fourth power, as is the case for annihilations to ZZ or 
W + W~ . For this reason, annihilations to ZHi, ZH 2 and W^H 1 " dominate over those to ZZ 
or W + W~ for neutralinos with a substantial bino fraction. 

As the expressions for neutralino annihilations to final states including Higgs bosons are 
very lengthy, we do not reproduce them here. For further details, we refer the reader to Ref. 
and references therein. 
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D. Neutralino Annihilations to 77 and 7Z 

Although neutralinos do not annihilate at tree level to 77 or 7Z, a large number of one-loop 
diagrams contribute to these processes [If. Although the cross sections for annihilations to 
these final states are considerably smaller than to heavy fermions, gauge or Higgs bosons, they 
are important because they produce mono-energetic gamma-ray lines, potentially providing a 
clear and distinctive signature in indirect detection experiments. Neutralino annihilations to 
77 and 7Z result in gamma-ray lines with energies of E 1 = m x o and E 1 = m x o — (m|/4m x o), 
respectively. 

Although the expressions for the cross sections to these final states are too long to reproduce 
here, we will briefly point out a few common features. In most cases, the largest contributions 
come from diagrams with a chargino-W 7111 loop. Each term in the amplitude from these processes 
is proportional to one of the higgsino fractions (liV^I 2 , liV^I 2 ), to the wino-fraction (|A" 12 | 2 ), or 
to some combination of these three {Ni^N* A , etc.). The bino content of the lightest neutralino 
does not contribute to this class of diagrams. Roughly speaking, we expect the annihilation 
cross section to gamma-ray lines to scale with the square of the higgsino fraction (or possibly 
the wino fraction) of the lightest neutralino. Even a pure-higgsino or pure-wino, however, has 
an annihilation cross section to 77 and jZ of around or less than 10 -28 cm 3 /s, which is less 
than 1% of the value needed to generate the observed thermal abundance. Lines are, therefore, 
always a small fraction of neutralino annihilation final states. 

In addition to diagrams with a chargino-H^ ± loop, diagrams involving a sfermion-fermion 
loop can generate the leading contribution to gamma-ray lines, in particular in the case of light 
sfermions. 

In Fig. |2l we plot the neutralino's annihilation cross section to gamma-ray lines (in the low 
velocity limit) as a function of higgsino composition. In models with heavy squarks (black points 
each represent models with rriq > 1 TeV), these cross sections are almost always dominated by 
W^-chargino loop diagrams, and therefore a strong correlation emerges between the higgsino 
fractions and the cross sections to 77 and jZ. This correlation is less strong among models 
with light squarks. These results, and the results shown in the other scatter plots included in 
this paper, were calculated using the DarkSUSY package @. 

E. Neutralino Annihilation Summary 

To summarize the content of this section, we will here briefly discuss under what circum- 
stances the various above mentioned neutralino annihilation modes are likely to dominate. To 
study this more concretely, we have plotted in figures El and 0] the fraction of neutralino annihi- 
lations which go to various final states as a function of the neutralino's bino-content. Each of 
the frames in these figures follow a contour of constant mass, such as those shown in figure ^ 

As a first case, consider a somewhat light neutralino, such that Higgs bosons, gauge bosons 
and top quarks are not kinematically accessible final states. Such a neutralino will annihilate 
almost completely to bottom quarks and tau leptons. If the sbottoms and staus are sufficiently 
heavy to not contribute, the leading annihilation channel will be through A exchange. The ratio 
of bottoms to taus produced through this channel is simply the ratio of the masses squared and 
the color factor, 3m 2 /m 2 ~ 10. Although interference terms from Z and A exchange diagrams 
can modify this ratio somewhat, annihilations to bottom quarks with a tau admixture of a 
few to ten percent is robustly predicted in the case of a light neutralino with much heavier 



8 



A 
V 



1e-28 



1e-30 



1e-32 



1e-34 



1e-36 




1e-28 



1e-04 0.001 0.01 0.1 

|N 13 | 2 



A 
V 



1e-28 



1e-30 



1e-32 



1e-34 



1e-36 





1e-28 



JO 



A 

a* 

V 



1e-30 - 



1e-32 - 



1e-34 



1e-36 



■ . " ■ i ■ "V 

■ ■ Si 

■-._f .»4C--»~."V 
' ■ " ■■v-in.:.-*' 



1e-10 



1e-10 1e-08 1e-06 1e-04 0.01 

|N 14 ! 2 



A 
> 



V 



1e-28 



1e-30 



1e-32 



1e-34 



1e-36 



■ ■ _ B l ■ ■ ■ ■ ■ ■ 

■ +■ ■% _■ . ." ■■ ■ ■ ■ 



1e-04 



0.001 



0.01 

|N 13 | 2 +|N 14 ' 2 



1e-38 



1e-04 



0.001 



0.01 

|N 13 | 2 +|N 14 | 2 




FIG. 2: The neutralino annihilation cross section to 77 and 7Z final states as a function of its higgsino 
content, | A^i3 1 2 , |A r i4| 2 and | iVi3 1 2 + | A r i4 1 2 . Dark (black) points represent models with heavy squarks 
(rriq > 1 TeV), while lighter (red) points have no such constraint. Each point shown represents a set of 
parameters within the MSSM that is not in violation of direct collider constraints and that generates 
a thermal relic density within the la range measured by WMAP (0.119 > tt x o > 0.0787) j5(] . The 
GUT relationship between the gaugino masses has been adopted. 
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sfermions. 

If the sfermions are not so heavy, however, and sfermion exchange diagrams to fermion pairs 
are substantial, then the ratio of annihilations producing bottoms and taus can be much larger 
or smaller. This depends on the sbottom and stau masses and on competition between sfermion 
exchange diagram and various interference terms. For example, in figures E3 and 01 changing 
the sfermion masses from very heavy (10 TeV) to 300 GeV 1 has the effect of slightly increasing 
the fraction of annihilations which produce taus. 

We next consider a neutralino heavier than the gauge bosons. In this case, since annihilations 
to W + W~ and ZZ have cross sections which approximately scale with ( | iV"i3 1 2 — |iVi4| 2 ) 2 , 
neutralinos annihilate primarily to fermions unless their higgsino fraction is rather large. 

If the lightest neutralino is heavy enough to annihilate to Higgs boson pairs, or to a gauge 
boson along with a Higgs boson, these modes are often important. As said before, annihilations 
to ZHi and ZH 2 are in cases proportional to Nf 3 rather than to the forth power, as is the 
case for annihilations to ZZ or W + W~ . For bino-like or mixed bino-higgsino neutralinos, 
these annihilations (especially those to ZHi, ZH 2 and W ± H T ) dominate over those to ZZ 
or W + W~ . Whether these modes occur more often than those to fermion pairs depends on 
features such as the value of tan/3, the sfermion and A masses, and the precise composition 
of the neutralino. For heavy neutralinos, annihilations to top quarks can also be significant, 
especially in models with low or moderate tan (3. 

In addition to the annihilation fractions, in figures |3] and H] we plot the thermal relic abun- 
dance of the lightest neutralino, and compare this with the quantity of cold dark matter mea- 
sured by WMAP @ (0.119 > Vt CDM > 0,0787 at the 2a confidence level). This quantity was 
calculated using the DarkSUSY package 0. 



1 Not all sfermion masses are set to 300 GeV and 1000 GeV in figures U3 and 0] once splitting has been taken into 
consideration. Although we have minimized the mass splitting by adopting A t .b. r soft terms much smaller 
than the sfermion masses, off diagonal terms in the sfermion mass matrices are proportional to /jm// tan/3 
for up-type fermions and /irrif tan for down type fermions, and can be large for bino-like neutralinos. 
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FIG. 3: The dominant neutralino annihilation modes for m^o =100 and 500 GeV, as a function of its 
bino content, tan/3 has been set to 10. The GUT relationship between Mi and Mi has been used, 
and in each frame tan/3 = 10 and 771^=300 GeV. In the left frames, 300 GeV and 1 TeV sleptons and 
squarks have been used, while in the right frame they have been made very heavy for illustration. Also 
shown above each frame is the relic abundance of the lightest neutralino compared to the quantity 
measured by WMAP. See the text for more details. 
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FIG. 4: The same as in figure EJ but with tan/3 = 50. See the text for more details. 
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IV. NEUTRALINO-NUCLEI ELASTIC SCATTERING 



Neutralinos scatter elastically with nuclei through both scalar and axial-vector couplings. 
Before moving on to direct and indirect detection methods, we review here the relevant elastic 
scattering cross sections for a neutralino. 



A. Scalar (Spin-Independent) Scattering 



A neutralino scattering elastically with a nucleus through scalar interactions does so with a 
cross section of: 



&si 



Am 2 x0 m 2 A 



: [Zf p + {A-Z)f n 



11) 



7r(m x o + vtiA) 2 1 

where rriA, Z and A are the mass, atomic number and atomic mass of the target nucleus. f p 
and /„ are the neutralino couplings to protons and neutrons, given by: 



fp,n — ^ ] /; 



m 



with f$ 



0.020, f$ 



q=u,d,s 

0.026, 



m, 



p,n 



0-118, f$ 



q=c,b,t 

0.014, /§> 



(12) 



0.036 and « 0.118. 

These quantities are subject to uncertainties in the relevant nuclear physics measurements (for 
more information, see Ref. flol|). a q are the neutralino-quark couplings. The first term in 
this expression corresponds to interactions with the quarks in the nucleus, either through t- 
channel CP-even Higgs exchange, or s-channel squark exchange. The second term corresponds 
to interactions with the gluons in the nucleus through a heavy quark/squark loop diagram. ffQ 
is given by 1 - fg - fg - fg 



0.84, and analogously, fi n ^ -~ 



1 TG 



0.83. To account for finite 
momentum transfer, the calculation should also include the appropriate form factor. 

The neutralino-quark coupling, in which all of the SUSY model-dependent information is 
contained, is given by [ill ]: 
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where throughout i = 1 for up- type quarks and i = 2 for down type quarks. mij,m 2 i denote 
elements of the appropriate 2x2 squark mass matrix and r] is the matrix which diagonalizes 
that matrix, T 3i and denote hypercharge, isospin and electric charge of the quarks. 
For scattering off of up-type quarks: 

5\ = Ni 3 , 82 = Nu, B — sin/3, C = sina, D = cosa, (15) 

whereas for down-type quarks: 

Si = ./V14, 82 = — iVi3, 5 = cos/3, C = cosa, D = — sina. (16) 

Here, a is the Higgs mixing angle. 

The first two terms of Eq. EH correspond to interactions through the exchange of a squark, 
while the final term is generated through Higgs exchange. As in the treatment of the anni- 
hilation cross section to fermions, the cross section resulting from Higgs exchange terms is 
proportional to the higgsino fraction (or possibly the difference of the two higgsino fractions) 
and bino fraction, thus favoring a mixed neutralino. Squark exchange diagrams, on the other 
hand, contribute for either a mixed neutralino or a pure bino. 

Taking a closer look at the contribution from Higgs exchange, note that if ^ 
mff 2) then cos a ~ 1 and sin a « 0. In this limit, the above expressions simplify con- 
siderably, yielding a q /m q rs gig 2 NnN* 4 /4mw sin Prnfj for up-type quarks and a q /m q w 
—gig2NiiN* 3 /4mw cos Prrijj for down-type quarks. The cross section therefore scales with 
terms such as |-/Vn| 2 |A^i4| 2 /m^ and | iVn | 2 1 iV"i3 1 2 tan 2 f3/m A Hi . The contributions from squark 
exchange can be similar in importance and, in the case of nearly diagonal squark mass ma- 
trices, scale as \Nn\ 2 \Ni^\ 2 /m~ for up-type quarks and |iVii| 2 |iVx3| 2 tan 2 /3/m~ for down- type 
quarks If any of these contributions are similar in magnitude, interference terms can also 
be important. 

In figures E] and El we plot the scalar elastic scattering cross section of the lightest neutralino 
as a function of its bino-content for various choices of tan /?, uia, fn x o and sfermion masses. 

B. Axial- Vector (Spin-Dependent) Scattering 

Although current direct detection experiments are primarily sensitive to spin-independent 
scattering, the elastic scattering of neutralinos with nuclei through axial-vector couplings is 
relevant to the capture rate of neutralinos in the Sun, and the corresponding indirect detection 
rate in neutrino telescopes. With this in mind, we will briefly summarize the most important 
aspects of this cross section here. 

Axial-vector elastic scattering can be induced by the t-channel exchange of a Z boson and 
by the s-channel exchange of a squark. The latter of these processes always generates a small 
cross section, (<Jsd ~ 10 -6 pb or less with a nucleon), however, which is well below the reach 
of planned direct detection experiments, and well below the magnitude needed to generate an 
observable flux of neutrinos from the Sun. Elastic scattering through Z exchange, however, can 
be much more significant. 

The axial-vector elastic scattering cross section of a neutralino with a nucleus through a Z 
boson simply follows from the Xi ~~ Xi ~ Z coupling and is proportional to cjsd oc (|Ai 3 | 2 — 
I iV"!4 1 2 ) 2 . This can lead to a cross section not far below ~ 10~ 2 pb in extreme cases. As we will 
discuss in section IV Bl this cross section can potentially generate an observable flux of neutrinos 
from the Sun if |A 13 | 2 - |A 14 | 2 is a few percent or larger. 
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m x o=100, m A =300, 300 GeV sfermions 



m x »=100, m A =300, heavy sfermions 





FIG. 5: The scalar (spin-independent) neutralino-nucleon elastic scattering cross section for m x o =100 
and 500 GeV, as a function of its bino content. The GUT relationship between M\ and M2 has been 
used, and in each frame tan/3 = 10, 50 (as labeled) and m,4=300 GeV. In the left frames, 300 GeV and 
1 TeV sleptons and squarks have been used, while in the right frame they have been made very heavy 
for illustration (10 TeV). Current bound on this quantity from the CDMS experiment is 2 x 10~ 7 pb 
or 7 x 10~ 7 pb for a 100 GeV or 500 GeV neutralino, respectively ^3). Also shown above each frame 
is the relic abundance of the lightest neutralino compared to the quantity measured by WMAP. See 
the text for more details. 
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m x o=100, m A =1200, 300 GeV sfermions 



m x <>=100, m A =1200, heavy sfermions 
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V. DEDUCING DARK MATTER PROPERTIES FROM ASTROPHYSICAL OB- 
SERVABLES 



The first step toward identifying the features of the supersymmetric model represented in 
nature (the lightest neutralino's composition, tan/5, etc.) is to determine the relationships 
between observable astrophysical quantities and the neutralino's phenomenological character- 
istics. In our study, we will consider the following astrophysical observables: 

• The rate in direct dark matter detection experiments 

• The flux of high-energy neutrinos from the Sun, generated through neutralino annihila- 
tions 

• The shape of the gamma-ray spectrum from neutralino annihilations (including lines) 

• The annihilation rate of neutralinos in the local halo, as inferred by the cosmic positron 
flux 

These can, in turn, potentially be used to infer the following phenomenological characteristics 
of the lightest neutralino: 

• The spin-independent elastic scattering cross sections of the lightest neutralino with nu- 
cleons 

• The spin-dependent elastic scattering cross sections of the lightest neutralino with nucle- 
ons 

• The relative fraction of neutralino annihilations which result in 77 or 7Z final states (at 
low velocities) 

• The magnitude of the lightest neutralino's annihilation cross section (at low velocities) 

Using this set of astrophysical observables, it may be possible to infer some of the properties of 
the underlying supersymmetric model. In this section, we will consider how these phenomeno- 
logical quantities may be measured in near future astrophysical observations, and what these 
measurements may be able to tell us about the nature of supersymmetry. 

A. The Neutralino's Elastic Scattering Cross Section 

There are several experiments currently in operation that hope to detect neutralino dark 
matter particles by observing them elastically scattering off of a detector. This class of tech- 
niques, known as direct detection, has become very exciting in recent years, reaching the level 
of sensitivity needed to explore regions of the supersymmetric parameter space. Furthermore, 
a large fraction of the remaining supersymmetric parameter space is within the projected reach 
of direct detection experiments planned for the near future. 

Currently, the strongest constraint on the neutralino's spin-independent elastic scattering 
cross section come from the CDMS experiment 0|. For a 50-100 GeV neutralino, CDMS 
has excluded <j x n > 2 x 10~ 7 pb, assuming a standard local dark matter density (p x o pa 0.3 
GeV/cm 3 ) [14]. For a 1 TeV neutralino, this constraint is about an order of magnitude weaker. 
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FIG. 7: The relationship between the quantity |iVn| 2 |iVi3| 2 tan 2 (5/m\ and the spin-independent 
neutralino-nucleon elastic scattering cross section. For models with a large cross section (> 10 
pb), the correlation is quite strong. Each point shown represents a set of parameters within the 
MSSM that is not in violation of direct collider constraints and that generates a thermal relic density 
within the 2a range measured by WMAP (0.119 > O x o > 0. 0787)0. The GUT relationship between 
the gaugino masses has been adopted. 



Constraints from the Edelweiss Zeplin jig] and WARP [l7j experiments are not very much 
weaker than those of CDMS. 

The spin-independent neutralino-nucleon elastic scattering cross section is a function of many 
super symmetric parameters, including the composition of the lightest neutralino, the masses 
and mixings of the exchanged Higgs bosons and squarks, and tan f3. This can make it difficult 
to extract much information about the supersymmetric model from a measurement of this 
quantity. If the cross section is measured to be quite large, however, it would indicate that 
this process is likely dominated by the exchange of the heavy CP-even Higgs boson, Hi, with 
a large value of tan/5 and a non-negligible higgsino fraction. In this limit, the cross section is 
roughly proportional to: a X N oc |-/Vii| 2 |Aq 3 | 2 tan 2 (3/m\ (see section ITV A|) . In Fig. we plot 
the relationship between these quantities, and demonstrate that in models with a large cross 
section (~ 10 -7 pb or larger) that the correlation is very strong. 

If a positive detection is made by CDMS or another direct detection experiment, such a 
measurement could be used to constrain both the neutralino's elastic scattering cross section 
and its mass. With an accumulation of ~ 100 events, the neutralino's mass could be measured 
to ~ 25% accuracy, and the cross section even more accurately (up to uncertainties in the local 
halo profile) [1S|] . With fewer events, these parameters will be considerably less constrained. 
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B. The Flux of High-Energy Neutrinos From The Sun 



Neutralinos traveling through the solar system can elastically scatter off of massive bodies 
such as the Sun, potentially becoming gravitationally bound. If this occurs sufficiently fre- 
quently, the neutralino annihilation rate in the Sun's core may be large enough to generate an 
observable flux of high-energy neutrinos 3 . 



Unlike rates from other methods of indirect detection, the neutrino flux from the Sun depends 
primarily on the neutralino's elastic scattering cross section with nuclei in the Sun rather than 
on its annihilation cross section. The annihilation cross section drops out of the calculation 
in the case that equilibrium is reached between the solar capture and annihilation rates. A 
measurement of the neutrino flux from the Sun, therefore, can be used to infer the elastic 
scattering cross section of the lightest neutralino, much as direct detection experiments do. 
Unlike the event rates produced in CDMS, Zeplin or Edelweiss, however, the neutralino capture 
rate in the Sun is the result of both axial- vector (spin-dependent) and scalar (spin-independent) 
interactions. Although direct detection experiments constrain the spin-independent elastic 
scattering of a neutralino, spin-dependent scattering is not yet constrained at a practical level. 
Furthermore, since it will be a number of years before next generation neutrino telescopes 
IceCube |2o[ and Antares |2l[ are fully operational, the limits on spin-independent scattering 
are likely to becomes stronger by at least an order of magnitude and probably more by the 
time the neutrino flux from the Sun will be measured. In models in which direct detection 
experiments will not discover neutralino dark matter in the near future, neutralinos will be 
potentially observable at next generation neutrino telescopes only if they are captured in the 
Sun primarily though spin-dependent elastic scattering. By combining information provided 
by direct detection experiments with the rates from the Sun observed by neutrino telescopes, 
it may be possible to measure the neutralino-proton spin-dependent elastic scattering cross 
section. 2 

The spin dependent neutralino-proton elastic scattering cross section is generated through 
two classes of diagrams: t-channel Z exchange and s-channel squark exchange. The contribution 
from squark exchange, however, is invariably small (typically 10~ 6 pb or less) compared to the 
value needed to generate an observable flux of neutrinos from the Sun. For any model that 
is potentially observable in existing or planned neutrino telescopes, neutralino spin-dependent 
elastic scattering will, therefore, be dominated by Z exchange. This cross section, following 
from the x\ ~ Xi ~ Z coupling, is proportional to |iV 13 | 2 — |iV 14 | 2 . In Fig. |SJ we plot the 
relationship between the neutralino-proton spin-dependent elastic scattering cross section and 
the quantity |iVi 3 | 2 — |Ai 4 | 2 . 

In figure 03 we plot the event rate in a kilometer-scale neutrino telescope as a function 
of the spin-dependent neutralino-proton elastic scattering cross section and as a function of 
1-^13 1 2 — | iV"i4 1 2 . We can see from this figure that a measurement of the high-energy neutrino flux 
from the Sun can be used as a rough measurement of the quantity |-/V 13 | 2 — |iV 14 | 2 . Furthermore, 
unless the higgsino fraction of the lightest neutralino is very large, the quantity |iV 13 | 2 — |A" 14 | 2 
translates fairly accurately to the neutralino's higgsino fraction, | iV 13 | 2 + |iV 14 | 2 . This is demon- 



2 It is also possible that future direct detection experiments will be sensitive to spin-dependent neutralino 
elastic scattering. The COUPP collaboration, for example, is developing bubble chamber detector technology 
that could potentially reach the level of sensitivity needed to probe parameter space of the MSSM through 
spin-dependent neutralino scattering [2^ . 
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FIG. 8: The relationship between the neutralino-proton spin-dependent (axial-vector) elastic scattering 
cross section and the quantity | ^Vi3 1 2 — | iV"i4 1 2 . For models along the diagonal line extending from the 
lower left to upper right portions of the figure, the process is dominated by Z exchange. Those points 
scattered from this line obtain their spin-dependent elastic scattering cross section significantly from 
squark exchange. Each point shown represents a set of parameters within the MSSM that is not in 
violation of direct collider constraints and that generates a thermal relic density within the 2a range 
measured by WMAP (0.119 > O x o > 0.0787)0]. The GUT relationship between the gaugino masses 
has been adopted. 

strated in figure in this way, neutrino telescopes may provide a rough measurement of the 
lightest neutralino's higgsino fraction 23]. By studying the spectrum of neutrinos, information 
regarding the lightest neutralino's mass may also be inferred |24| . 

Measurements of the neutrino flux from the Sun will be limited by the background of atmo- 
spheric neutrinos. Atmospheric neutrinos generate roughly ~ 500 muons per square kilometer 
per year above 50 GeV in the angular window corresponding to the Sun. Over ten years, a 3a 
detection would therefore require 10 x R v > 3 x y/10 x 500, or an annual rate in IceCube or 
KM3 of ~ 20 events. If a muon energy threshold of 100 GeV were imposed, a rate of 8 events per 
year would constitute a 3-cr detection. To be observable in a kilometer-scale neutrino telescope 
such as IceCube (or a kilometer-scale version of Antares, such as KM3), an event rate on the 
order of ~10 events per square kilometer per year will be required. If such a rate is observed, 
it would indicate a neutralino with a substantial higgsino component (~1% or greater). The 
absence of any signal of this magnitude could be used to constrain |Aq 3 | 2 — |A" 14 | 2 < 0.05. 

C. The Shape of the Gamma Ray Spectrum 

Gamma-rays can be produced in dark matter annihilations through a variety of annihilation 
modes. If annihilations occur sufficiently frequently in nearby regions, these gamma-rays may 
be observable by operating or upcoming Atmospheric Cerenkov Telescopes such as MAGIC [25], 
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FIG. 9: The rate of neutrinos detected in a kilometer scale neutrino telescope, such as IceCube, from 
neutralino annihilations in the Sun as a function of the neutralino's spin-dependent elastic scattering 
cross section (upper), and as a function of the quantity | iVi3 1 2 — | 1 2 (lower). While the current 



CDMS constraint on the spin-independent elastic scattering cross section 13] has been applied in the 
left frames, a constraint 100 times more stringer has been applied in the right frames in anticipation 
of increased sensitivity from direct detection experiments in the coming few years. Each point shown 
represents a set of parameters within the MSSM that is not in violation of direct collider constraints 
and that generates a thermal relic density within the 2a range measured by WMAP (0.119 > £l x o > 
0.0787)0]. The GUT relationship between the gaugino masses has been adopted. Each model shown 
contains a LSP heavier than 100 GeV, in order to avoid strong suppression from the muon energy 
threshold of 50 GeV that has been used. 



HESS |26| or VERITAS j2?|, or by satellite-based experiments such as GLAST j28|. Amon g th e 
most promising sources of dark matter annihilation radiation are the center of our galaxy [29], 
dwarf spheroidal galaxies (such as Draco, Sagitarius, or Canis Major) |3Cj |. external galaxies 
|3l| . galactic intermediate mass black holes |32(, and small scale dark matter substructure j33[. 

If one or more of these sources is identified, the spectrum of gamma-rays could (in principle) 
be used to learn about the modes through which neutralinos annihilate. In figure ^2 we plot 
the gamma-ray spectrum, per neutralino annihilation, for several of the annihilation modes 
likely to contribute significantly. From this figure, it is clear that most of annihilation modes 
result in a similar spectrum. The primary exception to this is the spectrum from annihilations 
to tau lepton pairs, which has a considerably harder slope than the other modes. The spectrum 
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|N 13 | 2 - |N 14 | 2 

FIG. 10: The higgsino fraction of the lightest neutralino (I-/V13I 2 + | A^i^ | 2 ) verses the quantity | iVi3 1 2 — 
|iVi4| 2 . For neutralinos with less than roughly 10% higgsino fraction, a measurement of | iVi3 1 2 — | -ZV14 1 2 
can act as a measurement of the neutralino's higgsino fraction. Each point shown represents a set of 
parameters within the MSSM that is not in violation of direct collider constraints and that generates 
a thermal relic density within the 2a range measured by WMAP (0.119 > U x o > 0.0787) ja]. The 
GUT relationship between the gaugino masses has been adopted. 

is also slightly harder for annihilations to W^H^. 

A gamma-ray telescope, at best, might hope to identify the fraction of neutralino annihilation 
which go to the distinctive t + t~ mode. The ability to make such a measurement, of course, 
depends on the brightness of the annihilation source (or sources). The number of gamma- 
rays from dark matter annihilations to be seen by GLAST is constrained by the lack of such 
detections by EGRET. Although GLAST will be considerably more sensitive than EGRET, 
it is unlikely that it will detect more than ~10-20 gamma-rays from a source of neutralino 
annihilation. If multiple such sources are discovered by GLAST, this number could be larger. 

In figure we show how difficult it is to identify a r + r~ component in the gamma-ray 
spectrum from neutralino annihilations. This is done by comparing the number of gamma-rays 
observed at low energies (1-15 GeV) to those observed at higher energies (above 15 GeV). For 
as neutralino annihilating purely to b quarks, the ratio of low energy to high energy gamma-rays 
is ~ 7.2, while for annihilations to tau pairs yields a much smaller ratio (~ 0.3). Many events 
will be required to accurately measure this quantity, unfortunately. The outer and inner error 
bars shown correspond to a total number of 15 and 100 events, respectively. With 15 events, 
a neutralino annihilating only to bb is indistinguishable from a WIMP annihilating most of 
the time to t + t~. The situation is not all that much better with 100 events. It appears that 
gamma-ray measurements will not be able to tell us much about the annihilation modes of a 
neutralino unless it annihilates almost entirely to taus. 

In addition to the continuum emission, monoenergetic gamma-ray lines can be generated 
through loop diagrams in which neutralinos annihilate directly to final states including photons, 
X1X1 ~~ y 77 an d X1X1 ~ y 1% [8]. The cross sections to these final states are only a small fraction 
of the total neutralino annihilation cross section, however, making lines more difficult to observe. 
If the dark matter density is large in the galactic center region, (ie. a cusped j34| or adiabatically 
compressed [35] profile, or a density spike resulting from the adiabatic accretion of dark matter 
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FIG. 11: The gamma ray spectrum per neutralino annihilation for various modes and two choices 
of the neutralino mass. With the exception of the t + t~ mode (and to a lesser extent, W^H^), 
the spectrum produced is very similar. The harder spectrum produced through taus is potentially 
distinguishable from other modes, however. Masses of rn^ Hl H ± ~ 300 GeV and m# 2 ~ 115 GeV 
were used. 
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FIG. 12: The ratio of gamma-rays above and below 15 GeV, as a function of the annihilation fraction 
to t + t~ (with the remaining annihilations to bb), for a 500 GeV neutralino. The large and small error 
bars are for a total number of events of 15 and 100 in GLAST, respectively. A very large number of 
events would be needed before the fraction of neutralinos to t + t~ could be measured. 
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onto the central supermassive black hole [36]) it is possible that such lines could be detected. 
Observations of this region by HESS j3?| and MAGIC 38|| have revealed a background from this 
region which will make gamma-ray lines difficult to detect, however [39]. In light of this, other 
regions (such as dwarf spheriodal galaxies |30() appear to be more likely sources of observable 
neutralino annihilation lines. 

The observation of gamma-ray lines from neutralino annihilations would provide both a 
measurement of the neutralino's mass, and information regarding the neutralino's composition. 
Although the processes which lead to 77 and 7Z final states include a wide variety of loop 
diagrams (involving sfermions, charginos, neutralino and Higgs bosons), making it difficult to 
determine the responsible supersymmetric parameters from so few observables, we have shown 
in Fig. 12 that neutralinos with a substantial higgsino (or wino) component generate brighter 
lines than highly bino-like neutralinos. 



D. The Cosmic Positron Flux 

If gamma-rays are detected from neutralino annihilations, they will not tell much about the 
magnitude of the responsible cross section. This is because little is known about the density of 
dark matter in the inner region of halo profiles, which is critically important in determining the 
resulting annihilation rate. The annihilation rate anticipated in the central region of our galaxy, 
for example, can vary by at least ~8 orders of magnitude depending on whether a profile with a 
flat core is present, or instead an adiabatically compressed cusp j35[, or a density spike generated 
through the accretion of dark matter onto the galaxy's supermassive black hole 36] . Given this 
uncertainty, a measurement of the annihilation rate in the region of the galactic center will tell 
us little about the neutralino's annihilation cross section. The uncertainties associated with 
the halo profiles of other objects (dwarf spheroidal galaxies, for example), although less severe, 
are still sufficiently large to make cross section measurements impractical. 

If the annihilation rate could be measured in a region without such enormous uncertainties in 
the halo profile, however, perhaps a determination of the neutralino's annihilation cross section 
could be made. The cosmic positron flux is, perhaps, the most useful tool for making such a 
determination. 

Unlike with gamma-ray measurements, cosmic positron (as well as anti-proton) measure- 
ments potentially observe the annihilation products of dark matter produced over large volumes 
of space (several kpc 3 ). Such a measurement, therefore, can be used to determine the product of 
the neutralino's annihilation cross section and the neutralino density squared, averaged over the 
sampled volume. The volume sampled roughly corresponds to the distance a typical positron 
travels from its point of origin before losing the majority of its energy; a few kiloparsecs for 
positrons in the relevant energy range. As a result of this limited range, only the dark matter 
distribution in the local halo is relevant to the observed cosmic positron flux. Assuming there 
are no very large and unknown clumps of dark matter in the surrounding kiloparsecs (which, 
although not impossible, is very unlikely 46]), a measurement of the cosmic positron spectrum 
could be used to infer the neutralino's annihilation cross section (in the low velocity limit) with 
only a relatively modest degree of astrophysical uncertainty. 

Future measurements of the cosmic positron spectrum will potentially reach the level of 
sensitivity needed to detect a contribution from neutralino annihilations [47[ . In Figs. and 
IT4^ we show the ability of the AMS-02 experiment to detect the presence of positrons generated 
in neutralino annihilations. 
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FIG. 13: The positron spectrum for a 100 GeV neutralino annihilating to bb (left) and to W + W~, 
ZZ (right). Also shown as error bars is the projected ability of AMS-02 to measure this spectrum. 
An annihilation cross section of 3 x 10 -26 cm 3 /s, and a boost factor of 5 have been used. 
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FIG. 14: The x 2 over the expected background of a detection of cosmic positrons generated in neu- 
tralino annihilations, as projected for AMS-02. Along the x-axis is the neutralino annihilation cross 
section multiplied by the boost factor resulting from local inhomogeneities in the dark matter dis- 
tribution. For somewhat light neutralinos, < av > x BF > 2 x 10 -26 cm 3 /s could be detected by 
AMS-02 El. 



VI. INPUT FROM COLLIDERS 



If TeV-scale supersymmetry exists in nature, it will very likely be discovered at the Large 
Hadron Collider (LHC), or perhaps even earlier at the Tevatron. Although the work presented 
here is not a collider study, we will briefly discuss here some of the measurements that could 
be made at the LHC or Tevatron that would be complementary to the information ascertained 
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through astrophysical experiments. For further discussions regarding neutralino dark matter 
at collider experiments see, for example, Refs. H, H, 4t| 5(], 51, 12, 53l . 



Of the searches for new physics to be conducted at the LHC, those for squark and gluinos 
will be among the most easy to perform. Squarks and gluinos are each produced with very large 
cross sections jH^I, and their subsequent decays result in distinctive jets plus missing energy 
signatures (assuming R-parity is conserved). After about one month of running (at the first 
year design luminosity), the LHC should be able to discover squarks or gluinos if either are 
lighter than about 1.5 TeV. Of course, this is dependent upon understanding the behavior of 
the experiment, which could very well take a year or so to achieve. Ultimately, the LHC will 
be sensitive to squarks and gluinos with masses up to ~ 3 TeV j55[. 

By studying the decays of squarks and/or gluinos, it will also be possible to discover 
others superpartners at the LHC. For example, in many models, decays of the variety, 

q -> x\q -> 

q — > xV- + l q, provide a clean signal of supersymmetry in the form of 
l + l~ + jets + missing E T . By studying the kinematics of these decays, the quantities nig, 
m x o , mi and m x o can each be reconstructed jiH HtJ , that is if the sleptons are sufficiently 
light for this decay chain to take place. More generally speaking, the LHC is in most models 
likely to measure the mass of the lightest neutralino to roughly 10% accuracy, and may also be 
able to determine the masses of one or more of the other neutralinos, and any light sleptons 
[55l | . Charginos are much more difficult to study at the LHC. 

Production of heavy, neutral Higgs bosons (A, Hi), can also be studied at the LHC. In 
particular, in models with large tan ft, heavy Higgs bosons have enhanced couplings to down- 
type fermions, thus leading to potentially observable di-tau final states. If enough of these 
events are observed, the masses of the heavy Higgs bosons could be potentially reconstructed, 
and tan ft measured H, 5^ . 

Prospects for the discovery of supersymmetry at the Tevatron, although not nearly as strong 
as at the LHC, are also very exciting. The most likely discovery channel at the Tevatron is 
probably through clean tri-lepton plus missing energy events originating from the production 
of a chargino and a heavy neutralino, followed by a decay of the form, X X2 ~ * vl^l^F — » 
vxV> ± l + l~x\ Only models with rather light gaugino masses (neutralinos and charginos) 
can be discovered in this way, however. For some of the recent results from supersymmetry 
searches at the Tevatron, see Ref. [Hl|. 

Determining the structure of an underlying supersymmetric theory from collider observables 
is by no means a simple task. Considerable effort has been recently been put into the "in- 
verse problem" jHsl IHH . I62I l63l . of extracting supersymmetric parameters from LHC data. 
Computational tools such as Sfitter |63[ and Fittino (64| have been developed for this purpose. 

We are not going to attempt a detailed collider analysis in this paper. Instead, we will simply 
adopt the position that, at a future time in which astrophysical dark matter experiments may 
be measuring the neutralino's phenomenological characteristics relevant to direct and indirect 
detection, collider experiments such as the LHC will be likely to have determined some of 
the basic features of the supersymmetric spectrum, such as the squark and gluino masses, the 
lightest neutralino's mass (and possibly one or more other neutralino masses) and possibly tan ft 
and rriA (in the case of large tan/3 and light ttla). In the more distant future, an International 
Linear Collider (ILC) could measure many more of these quantities, and at far grea ter precision 
[65] , revealing a great deal about supersymmetry and neutralino dark matter jl8| . 
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FIG. 15: An illustration of how supersymmetric parameters can be determined though astrophysical 
observations, aided by collider measurements. Solid lines represent methods that are more reliable or 
likely, whereas dotted lines are (qualitatively) less so. See text for further discussion. 



VII. PUTTING IT ALL TOGETHER: CONSTRAINING THE SUSY MODEL 

Thus far, we have discussed several ways of constraining supersymmetric parameters with 
astrophysical and collider measurements. We turn our attention now to combining these ob- 
servables in an effort to learn as much as possible about the characteristics of supersymmetry. 

In Fig. we illustrate how supersymmetric parameters, such as Mi and /i (or alternatively 
the mass and composition of the lightest neutralino), could be ascertained from astrophysical 
observations, aided by measurements at the Tevatron and/or LHC. In particular, we emphasize 
that information regarding the neutralino's composition is included in the elastic scattering 
rate, the neutrino flux from the Sun and in the brightness of gamma-ray lines. If we assist this 
process with measurements of the lightest neutralino's mass and squark masses at the LHC 
(as well as tan/3 and tua in some cases), we can combine these observables to hopefully arrive 
at a determination of supersymmetric parameters which are not well constrained by collider 
measurements alone. 

The complex interplay between these different observables is very difficult to discuss in a 
model independent fashion. To further study the collective reach of these experiments, we will 
consider a set of representative supersymmetric benchmark models. 
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VIII. BENCHMARK MODELS 



In an effort to evaluate how astrophysical measurements will be able to help constrain the 
nature of supersymmetry, in this section we will consider a series of benchmark models. We 
intend for these representative models to illustrate how various astrophysical observations could 
be used collectively to infer information about the composition of the lightest neutralino and 
other characteristics of supersymmetry. Although many sets of sup ersymmet r ic benchmark 
scenarios have been proposed elsewhere (see for example, Refs. |66l l67|). we have chosen to 
introduce our own here. The input parameters of these benchmark models are given in table |H 



Model 


M 2 




tan/3 


m f 


m A 


A t 


A b 






5a^ 




Notes 


LT1 
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55.7 
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3471. 
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2 
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817.9 
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10 
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3 
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10 
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TABLE I: The set of benchmark models we have adopted in our study. In addition to the input 
parameters of the model, the relic density and the amplitude of the contribution to the muon's magnetic 
moment a re g iven. Each point shown is in agreement with the constraint on the branching fraction 
B — > Xgj |68j as well as with all direct collider constraints. The labels LT, IM and HV denote models 
with light, intermediate mass and heavy neutralinos. 



In addition to the input parameters, in table we also show the relic density and the 
amplitude of the contribution to the muon's magnetic moment for each benchmark model, 
both calculated using DarkSUSY . For each model, the relic density falls within the 2a range 
measured by WMAP |5||. The value of the muon's magnetic moment, as measured using e + e~ 
data, is larger than the Standard Model prediction by 5a^ = (23.9 ± 7.2h a d-io ± 3.5ibi ± 6 cxp ) x 
10~ 10 , where the error bars denote theoretical uncertainties in the leading order hadronic and 
hadronic light-by- light contributions, as well as the experimental contribution |&|. Considering 
this data alone, our benchmark models LT2, IM1, IM2, HV1, HV2 and HV3 are each outside 
of the 2a range of this measurement. Given the unresolved issues involved in this measurement 
(such as the lower value found using r + r _ data) and the relativity low statistical significance 
of the measurement, we do not exclude any models on this criterion and leave it to the reader 
to account for this measurement as they choose. We have checked that each benchmark model 
is consistent with the measurement of the branching fraction B — > X s ^ [Hif • 
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IM5 


231.8 0.996 8 x 10" 5 0.003 0.0004 
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HV1 


413.8 0.996 6 x 10~ 5 0.003 0.0009 
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627.0 0.60 0.003 0.21 0.19 
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TABLE II: Phenomenological quantities for our benchmark models. From left to right are: the lightest 
neutralino's mass (in GeV), the bino, wino and higgsino fractions of the lightest neutralino, the spin- 
independent elastic scattering cross section of the neutralino with nucleons (in pb), the spin-dependent 
elastic scattering cross section of the neutralino with protons (in pb), the annihilation cross section of 
the lightest neutralino (in the low velocity limit, in cm 3 /s), the annihilation cross section to 77 and 
^fZ final states (in cm 3 /s), and the rate per square kilometer, per year in a neutrino telescope (with 
a 50 GeV muon energy threshold). 

In table EJ we show a number of phenomenological characteristics of the benchmark models, 
including the mass and composition of the lightest neutralino, as well as the most relevant 
elastic scattering and annihilation cross sections, and the rate predicted in neutrino telescopes. 
In the left frame of Fig. EE we plot the spin-independent elastic scattering cross sections of 
the lightest neutralino in these models compared to current and projected direct detection 
constraints. Although none of the benchmark models are excluded by current constraints, 
several of the models are within the projected reach of direct detection experiments in the near 
future. Experiments such as Super-CDMS (phase C), Zeplin-Max, EURECA, or XENON will 
be capable of testing all but the most heavy supersymmetric models (assuming that a neutralino 
makes up the dark matter of our universe). In the right frame of the figure, we compare the 
benchmark neutralinos' annihilation cross section to gamma-ray lines to the total annihilation 
cross section (in the low velocity limit). 

In table IIII1 we put forth a set of measurements which could plausibly be made by future 
experiments for each of our benchmark models. For the LHC, we focus on a few of the easiest 
to measure and most important quantities: m x o, m^, tan/3 and rriq. 3 The mass of the lightest 
neutralino can typically be measured to ~10% accuracy, as can squark masses if they are 
sufficiently light (in the case of very heavy squarks, we have estimated less precise measurements 



3 By niq, we mean the approximate masses of the first and second generation squarks. The stops and sbottoms, 
in particular, are sometimes significantly heavier and lighter than this value. 
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FIG. 16: In the left frame, we plot the spin-independent elastic scattering cross section for neutralinos 
in our benchmark models, compared to current and projected constraints. The solid line is the current 
constraint from the CDMS experiment The dashed line is the projected reach of CDMS (2007), 
whereas the dotted lines are (top to bottom) the projected reach of Super-CDMS Phase A, Zeplin-Max 
and Super-CDMS Phase C (see Ref. 7J|). In the right frame, we compare the benchmark neutralinos' 
annihilation cross section to gamma-ray lines (77 plus 7Z) to the total annihilation cross section. As 
dotted lines, we plot contours of constant o" 77+7 z/ctot ■ 



of this quantity), uia and tan/3 can potentially be measured if sufficiently light and large, 
respectively, by studying processes such as A — > t + t~. If this channel is observable, tan/3 is 
expected to be measured with a precision of approximately 15% It may be possible to 

determine rriA within a few GeV at the LHC. 

In some of these benchmark models, other supersymmetric particles may be measured at the 
LHC, especially in the lighter models. Often the masses of one or more additional neutralino 
and/or light slepton can be determined, for example. To explore this further would require a 
detailed collider study, which is beyond the scope of this paper. 

In addition to quantities measurable at the LHC, a number of astrophysical observables are 
considered in table 11111 Projecting the precision at which these quantities will be measured 
in the coming years is more difficult as they depend on a number of unknown astrophysical 
inputs. Measurements of the neutralino's elastic scattering cross section depend on the local 
dark matter density and, therefore, may be skewed if we happen to live in a dense clump of 
dark matter, or in a local void. Since these extreme deviations from the dynamically inferred 
average local density are quite unlikely we will proceed under the assumption that, with 
sufficient exposure, direct detection experiments will measure a X N to roughly within a factor 
of two. 

As we stated in section IV Bl the flux of high-energy neutrinos from the Sun will need to 
generate at least ~10 events per year if it is to be distinguished from the atmospheric neutrino 
background by experiments such as IceCube or KM3. This channel is even less subject to 
astrophysical uncertainties than direct detection because the capture rate of neutralinos in the 
Sun is averaged over the history of the solar system and, therefore, is robust to local fluctuations 
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TABLE III: A hypothetical set of measurements made by collider and astrophysical experiments for 
each of our benchmark models. Experiments considered include the Large Hadron Collider (LHC), 
direct detection experiments (CDMS, Zeplin, Edelweiss, and their planned extensions), neutrino tele- 
scopes (IceCube, KM3), antimatter detectors (AMS-02) and gamma-ray telescopes (GLAST, MAGIC, 
HESS, VERITAS, etc.). ± percentages and factors shown reflect the precision with which these mea- 
surements may plausibly be made. Entries reading as < X simply constrain the observable quantity. 
An entry of 'YES' indicates a measurement of a few GeV accuracy or better. Even if tan/3 and 
cannot be measured at the LHC (as denoted by 'NO' in the table), the lack of a detection can be used 
to constrain these quantities roughly as (tan/3/10)(200GeV/mA) ^ 1. 

in the dark matter distribution. 

Measurements of the cosmic anti-matter spectrum are subject to uncertainties in the galactic 
dark matter distribution and corresponding neutralino annihilation rate. Proceeding under the 
assumption that the positron boost factor is not surprisingly large (BF < 10), AMS-02 will 
either be capable of making a rough measurement of the neutralino's annihilation cross section 
(in the low velocity limit), or will place an upper limit on this quantity. 

It is impossible to know whether the inner cusps of the Milky Way and nearby dwarf sphe- 
riodal galaxies will be sufficiently concentrated to generate observable gamma-ray signals from 
annihilating neutralinos. If they do generate observable gamma-ray fluxes, then the brightness 
of 77 and 7Z lines, relative to the continuum emission, could be determined. If the brightness 
of lines is low compared to the continuum flux, this will be more difficult. 

We now will ask how much information regarding the supersymmetric parameters is gained 
by considering astrophysical measurements in addition to LHC data. We begin by considering 
the quantity fi. The parameter /i is difficult to constrain at the LHC in many supersymmetric 
scenarios. The spectrum and compositions of neutralinos and charginos could be studied to 
determine \i at a linear collider, but at the LHC it will likely go largely unconstrained. In con- 
trast, fi is very important to the phenomenology of neutralino dark matter, critically effecting 
the composition of the LSP and therefore its couplings and interactions. 



31 



In Figs. El and we plot the ability of the set of observations given in table IIHI to con- 
strain the parameter /x. For each benchmark model, we have randomly scanned over SUSY 
parameters 4 for models with observable quantities within the range of the measurements given 
in the table. For each frame (ie. for each benchmark model), the upper histogram reflects 
distribution of those models we found that satisfied the LHC measurements given in table ITTT| 
in addition to satisfying the relic density constraint. The lower histogram for each benchmark 
model are those models which also satisfy the astrophysical measurements given in table 11111 
The vertical dotted line in each frame denotes the actual value of /x for the given benchmark 
model. Alternatively, instead of the quantity /x, we could consider the ability of these mea- 
surements to constrain the lightest neutralino's higgsino fraction. Those results are shown in 
Figs. DH and E@ 

For many of the benchmark models, the inclusion of astrophysical data in the analysis allows 
for /x (and the higgsino fraction) to be significantly more constrained than by the LHC and relic 
abundance considerations alone. In a minority of the benchmarks, this additional information 
is not very useful, however. For model LT2, for example, the two histograms fall nearly on 
top of each other, indicating that astrophysical data does not allow for /x to be constrained by 
astrophysical data much further than with collider data (plus relic density) alone. Even in such 
a case, astrophysical data may be used to provide an important confirmation on the model 
assumptions that have been used. 

In five of our twelve benchmark models (IM1, IM2, IM5, HV1 and HV3), the quantities 
rriA and tan/3 are unlikely to be determined by the LHC. We have studied whether these 
quantities can be constrained by astrophysical observations and have plotted the histograms of 
the distribution of values of uia and tan (3 in Figs. ED and [221 (using the same methodology as 
Figs. fTTIOD|) . In benchmark models, IM1 and HV1, the lightest neutralino annihilates primarily 
through the A- resonance (the so-called "A-funnel" region of parameter space). In such models, 
rriA can be rather tightly determined if astrophysical data is taken into account. In a model 
such as HV3, on the other hand, little is gained by such information, tan/3 does not appear to 
be significantly constrained by astrophysical data in the models we have considered. 



We have scanned linearly over M2, fJ,, tua, tan/3, At, A^, and mj. M2, /x and mj have each been varied up to 
4 TeV, while vha and tan p have been varied up to 1 TeV and 60, respectively. A t and Af, have been varied 
up to 3 times mj. 

32 



LT1 





LT2 



i£_h 



-600 -400 -200 200 400 

H (GeV) 





600-600 -400 -200 200 400 

|i (GeV) 

LT4 ' 



600 



n 



1A 



-600 -400 -200 200 400 600-800 -600 -400 -200 200 400 600 800 

H (GeV) ji (GeV) 





o 



1000 



H (GeV) 



-1000 -500 500 1000 1500-1000 

H (GeV) 

FIG. 17: The ability of future astrophysical and collider (LHC) data to constrain the parameter \x in 
the first six of our benchmark models. In each frame, the upper histogram describes the distribution 
of models found in our random parameter scan which satisfy the LHC measurements given in table Hill 
in addition to the relic density constraint. The lower histogram for each benchmark model contains 
those models which also satisfy the astrophysical measurements given in table IIII1 The vertical dotted 
line in each frame denotes the actual value of fi for the given benchmark model. For models LT1, LT3, 
LT4, IM1 and IM2, astrophysical data allows for fi to be considerably more constrained. The model 
LT2 benefits far less from astrophysical data. The sign of [i is not expected to be determined by the 
techniques discussed here. 22 
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FIG. 18: The same as Fig. 1171 but for the last six of our benchmark models. For models IM3, IM4 
and HV3, astrophysical data allows for /i to be considerably more constrained. Models IM5, and HV1 
benefit less from astrophysical data. The value of \^\ is quite well constrained by LHC data alone (plus 
relic abundance) in model HV2. The sign of /i is not expected to be determined by the techniques 
discussed here. For model IM5, the lower curve has been multipled by 5 for illustration. 
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FIG. 19: The same as Fig. 1171 but showing the ability of future astrophysical and collider data to 
constrain the higgsino fraction of the lightest neutralino for the first six of our benchmark models. 
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FIG. 20: The same as Fig. I191 but for the last six of our benchmark models. For model IM5, the lower 
curve has been multipled by 5 for illustration. 
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FIG. 21: The ability of future astrophysical and collider (LHC) data to constrain the mass of the 
CP-odd Higgs boson, m^, in those of our benchmark models in which it is unlikely to be measured 
at the LHC. In each frame, the upper histogram describes the distribution of models found in our 
random parameter scan which satisfy the LHC measurements given in table II I II in addition to the 
relic density constraint. The lower histogram for each benchmark model contains those models which 
also satisfy the astrophysical measurements given in table UTTI The vertical dotted line in each frame 
denotes the actual value of txia for the given benchmark model. For models in the A-funnel region 
of paramter space (IM1 and HV1) astrophysical data allow for wia to be quite well constrained. The 
models IM2, IM5 and especially HV3 benefit fa^yless from astrophysical data. For model IM5, the 
lower curve has been multipled by 5 for illustration. 
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FIG. 22: The same as in Fig. I211 but reflecting the ability of future astrophysical and collider data 
to constrain tan/3 in those benchmark models that it is not determined at the LHC. In none of these 
models does astrophysical data substantially assist in determining the value of tan/3. For model IM5, 
the lower curve has been multipled by 5 for illustration. 
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IX. MODEL VARIATIONS AND OTHER CAVEATS 



Throughout this study, we have made a number of simplifying assumptions regarding the 
nature of supersymmetry. Our scans have been performed over the following set of parameters: 
M 2 , /i, rriA, tan/5, A t , Ab, and mj. M 2 , fi and mj have each been varied up to 4 TeV, while 
rriA and tan/3 have been varied up to 1 TeV and 60, respectively. A t and A^ have been varied 
up to 3 times mj. Positive and negative values of each parameter (except for tan/?) have been 
allowed. 

These seven parameters obviously do not span the entire 120 dimensional parameter space of 
the MSSM. We have reduced the number of free parameters considerably by adopting a single 
universal scalar mass (m ;) , and by fixing M\ and M3 to M 2 through the requirement that they 
unify at the GUT scale. Motivated by a desire to avoid flavor changing neutral currents and 
large electron and neutron electric dipole moments, we have taken all quantities to be diagonal 
in flavor space and have not considered the presence of any CP-violating phases. 

Of course we expect our results to change somewhat if any or all of these assumptions are 
relaxed. From the standpoint of the neutralino dark matter sector, the most dramatic variations 
occur when Mi is not related to M 2 by the GUT relationship. If M 2 is comparable to or smaller 
than Mi, for example, the lightest neutralino can be primarily wino-like or a mixed bino-wino. 
This is found in the case of Anomaly Mediated Supersymmetry Breaking (AMSB) models, for 
example, in which the gaugino masses are related by the ratios: Mi : M2 : M3 « 2.8 : 1 : 1.7, 
leading to the lightest neutralino being a nearly pure wino. The wino-like neutralinos found 
in AMSB models have a distinct phenomenology, annihilating very efficiently, including to 77 
and jZ final states (with nearly the largest cross section possible for a neutralino to these 
states) [73, E3 The large overall annihilation cross section leads to a thermal abundance well 



below the measured dark matter density, however, and thus AMSB scenarios require a non- 
thermal mechanism to generate the universe's dark matter [t3] - Assuming that such a process 
exists to generate the observed dark matter density in the form of wino-like neutralinos, this also 
leads to large cosmic positron and anti-proton fluxes F2L [73| , possibly capable of generatin g th e 
positron excess observed by the HEAT experiment |43j (and to some degree by AMS-01 |44j) 
without a large degree of local inhomogeneities [it. without a large boost factor) j4£j. The 
elastic scattering cross section of the wino-like neutralinos found in these models, in contrast, 
is somewhat smaller than is found in most other supersymmetric scenarios [t3 | . 

If the GUT relationship between Mi and M 2 is broken in the opposite way, such that Mi 
is less than M2/2, then the existence of a very light (10-50 GeV), bino-like neutralino becomes 
possible without violating the limits on charginos placed by LEP {m x ± > 104 GeV) |75j. If we 
look beyond the MSSM to models such as the NMSSM (the Next-to-Minimal Supersymmetric 
Standard Model), even lighter neutralinos are possible [76j. 

If, instead of adopting a single scale for the sfermion masses, we had selected values of 
rriQ, m?, m| m| and m? each independently (and independently for each family as well), 
the spectrum of squarks and sleptons could be considerably different from those considered in 
our study. In scenarios in which the scalar masses unify at a common scale, such as within 
the constrained MSSM (the CMSSM, or mSUGRA), itlq, m s and mg are typically similar is 
magnitude, where as and m s are sometimes smaller. This is a consequence of contributions 
to Wq, m\ and m| (but not to m| and m?) from the renormalization group running terms being 
proportional to the gluino mass. This leads roughly to m 2 L Pd m\ ~ m?Q — 8M|. Therefore, 
within the CMSSM, sleptons lighter than squarks are often predicted, particularly in the case 
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of sfermions lighter than or comparable in mass to Mi- The masses of the sleptons are not of 
particularly great importance for our study, however, as they do not effect the neutralino spin- 
independent or spin-dependent elastic scattering cross sections with nucleons, and typically 
contribute little to the neutralino annihilation cross section to lines. The most important effect 
of the sleptons is in the relic density calculation. If sleptons only slightly heavier than the 
lightest neutralino are observed at the LHC, for example, some of our conclusions may be 
modified. 

Of course there is another important assumption we have adopted throughout this study: 
that the spectrum of new particles being observed at colliders and in astrophysical experiments 
are supersymmetric particles. Even if a number of superpartners are discovered at the LHC, 
it is not at all clear that those particles will be identified conclusively as such. In particular, it 
has shown to be a challenge to distinguish supersymmetry from models with Universal Extra 
Dimensions (UED) at the LHC jz?]], especially if the Kaluza-Klein states appear with masses 
of ~ 1 TeV or above. Fortunately, the Kaluza-Klein dark matter found in these models [78| 
possesses features which cou ld p otentially be distinguished from neutralinos, such as a high 
neutrino rate from the Sun [79(, a hard positron spectrum [8Cj and fairly predictable direct 
detection rates [s^J - Furthermore, other scenarios, such as T-parity conserving little Higgs 
models may be challenging to distinguish from supersymmetry at the LHC, but might be 
assisted by dark matter observations. 



X. SUMMARY OF OUR RESULTS 

In this article, we have studied the ability of astrophysical dark matter experiments to 
constrain the parameters of supersymmetry. A number of observable quantities can be useful 
in this respect, including rates at direct detection experiments, the flux of high-energy neutrinos 
from the Sun, the brightness of gamma-ray lines and the flux of high-energy cosmic positrons. 

In exploring ways to use these observables to learn about the properties of supersymmetry, 
we have exploited a number of correlations: 

• In the case that the scalar (spin-independent) neutralino-nucleon elastic scattering cross 
section is relativity large (> pb), this process is dominated by the exchange of the 
heavy, CP-even Higgs boson, coupled to s and 6-quarks. In this case, the cross section 
scales as <j x n oc tan 2 /3\ Nu\ 2 \ Ni^\ 2 /m\. For smaller cross sections, contributions from 
light CP-even Higgs exchange and squark exchange can also play dominant roles, making 
the correlation less powerful, but still potentially useful. 

• As constraints on the spin-independent, neutrino- nucleon elastic scattering cross section 
become stronger, the only way that an observable flux of high-energy neutrinos from 
neutralino annihilations in the Sun could be generated is if neutralinos are captured in 
the Sun primarily through spin-dependent (axial-vector) couplings. In observable models, 
this is dominated by Z-exchange diagrams, leading to the strong correlation, oc 
(I ^13 1 2 — | 1 2 ) 2 . The rate in a neutrino telescope scales with this quantity as well (as 
long as the neutralino mass is well above the energy threshold of the experiment). 

• In most models, the neutralino annihilation cross sections to gamma-ray line producing 
final states (77, 7Z) are dominated by diagrams that include a chargino-PF 111 loop. In this 
class of diagrams, the cross section scales as some combination of the higgsino and wino 
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fractions of the neutralino. A measurement of the annihilation cross section to gamma- 
ray lines (or even the relative fraction of annihilations which produce lines) could be used 
to constrain the higgsino and wino components of the lightest neutralino. 

• The total neutralino annihilation cross section (in the low velocity limit) depends on a 
large number of supersymmetric parameters. Nevertheless, this quantity would have some 
discriminating power if it could be determined by measurements of the cosmic positron 
spectrum, which samples the dark matter annihilation rate in the local halo. 

Although the mass of the lightest neutralino will likely be measured to roughly 10% precision 
at the LHC, its composition will be considerably more difficult to constrain. To study the 
ability of astrophysical observations to assist in determining the composition of the lightest 
neutralino, we have considered a number of benchmark models. Asumming that the mass of 
the lightest neutralino and the squark masses (and in some cases tan f3 and tua) can be measured 
at the LHC, we performed parameter scans looking for models consistent with a given set of 
LHC observations, and consistent with the measured dark matter density. Comparing the 
distribution of the value of fi among these models (or alternatively, the higgsino fraction of 
the lightest neutralino) to the distribution after constraints from plausible future astrophysical 
observations were made, we find that in many scenarios, the astrophysical measurement can 
play an important role in determining the value of fi (or the higgsino fraction). 

In particular, we have found that: 

• In the majority of our benchmark models (LT1, LT3, LT4, IM2, IM3, IM4 and HV3), 
we found that the inclusion of information from astrophysical dark matter experiments 
(most importantly direct detection) allowed for a constraint on the magnitude of fi (or 
alternatively, the higgsino fraction of the lightest neutralino) that is considerably more 
stringent than can be obtained from collider experiments and relic density considerations 
alone. In addition to the important role of direct detection in constraining fi, we note 
that in those models which generate an observable neutrino flux from the Sun (IM2, 
IM4 and HV2) that \fi\ (and the higgsino fraction) are quite tightly constrained by this 
measurement. This is further enhanced by the very large annihilation cross sections 
neutralinos have to 77 and 7Z in these models. None of these measurements assist 
significantly in determining the sign of fi. 

• In a few of the benchmark models, astrophysical dark matter measurements provide 
far less information regarding the value of fi (and the higgsino fraction). For the case 
of model LT2, for example, nearly all of the supersymmetric models we found which 
were consistent with the projected measurements at the LHC were also found to predict 
similar elastic scattering cross sections and other astrophysically relevant quantities. For 
this reason, astrophysical measurements provide little discriminating power, although can 
still provide a valuable confirmation. Additionally, in the cases of models IM5 and HV1, 
little is learned about fi (or the higgsino fraction) from astrophysical measurements. 
In these cases, the LHC will likely not be capable of measuring tan/5 or m^, leading 
to a situation where a number of very different sets of parameters could lead to the 
measured relic abundance and astrophysical observables. This is further exasperated by 
the significant role of coannihilations in the freeze-out process of model IM5 and of the 
resonant annihilations through the CP-odd Higgs in model HV1. 
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• In two of our benchmark models (IM1 and HV1), the mass of the CP-odd Higgs could be 
quite constrained by astrophysical dark matter measurements. This is very interesting 
considering that in these models, A is too heavy and tan/3 is too small for heavy neu- 
tral MSSM Higgs bosons (A, Hi) to be detected at the LHC. In each of these models, 
neutralino annihilations in the freeze-out process is dominated by the resonant exchange 
of a CP-odd Higgs. When direct detection experiments determine that the lightest neu- 
tralino's elastic scattering cross section with nucleons is very small (7x 10 -11 and 2x 10~ 10 
pb, respectively), it demonstrates that these neutralinos are quite bino-like (a small hig- 
gsino fraction). This is further supported by the small annihilation cross section to 77 
and 7Z and the small rate in neutrino telescopes found in these models. A very bino-like 
neutralino requires an annihilation resonance (or possibly coannihilations) to avoid being 
overproduced in the early universe, thus it can be determined that rriA ~ 2m x o. In such 
models, tua can be constrained in this way to roughly ±100 GeV. 

• We have also studied the ability of astrophysical measurements to be used to constrain the 
quantity tan (3. We find that the prospects for this determination are far less optimistic. 

XI. CONCLUSIONS 

If low energy supersymmetry exists in nature, it will likely be discovered in the next few 
years at the LHC (or possibly at the Tevatron). Over roughly the same period of time, the 
prospects for direct and indirect searches for neutralino dark matter are also very encouraging. 
Each of these windows into the characteristics of supersymmetry can provide us with useful 
and complementary information. In this paper, we have studied the ability to constrain the 
parameters of supersymmetry, beyond what can be done at the LHC, using direct and indirect 
dark matter experiments. 

Although the presence of a neutralino LSP can in most cases be confirmed at the LHC, the 
composition and corresponding couplings of such a state will likely go unconstrained by such 
an experiment. Direct and indirect dark matter experiments, on the other hand, can probe the 
lightest neutralino's interactions with nucleons and with themselves (annihilations), potentially 
allowing for a determination of its composition (and the value of We also find that direct 
and indirect dark matter measurements, combined with relic abundance considerations, can 
in some cases (the A-funnel region) determine the mass of the CP-odd Higgs boson (A), even 
when beyond the reach of the LHC. 
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